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ieee! RODUC TION 


A. BACKGROUND 

Jet engines undergoing maintenance are tested in jet 
engine test cells prior to their operational use. A typical 
test cell is shown in Figure 1 and can be viewed as a "U" 
shaped concrete structure of rectangular cross-section. The 
cell is designed to provide clean distortion free airflow 
to the engine inlet, which is necessary for accurate 
performance evaluation. In addition, the cell attenuates 
noise and provides safety for test equipment and personnel. 

Engine exhaust is directed into an augmenter tube, then 
deflected vertically to be expelled into the atmosphere. 

The augmenter tube acts in combination with engine exhaust 
momentum to form a jet pump. In reality, the engine exhaust 
jet entrains secondary air (air not ingested by the engine) 
which is then carried through the augmenter tube by exhaust 
jet momentum. Frictional losses associated with the mass 

of air moving through the test cell are thus overcome by the 
jet pump's energy addition. 

The quantity of secondary air entrained is crucial to 
proper engine testing. Too little entrainment allows exhaust 
gas recirculation to the engine inlet. Too much entrainment 
causes excessive static pressure gradients to exist between 
engine inlet and exhaust planes. Static pressure should be 


nearly identical in both planes for proper testing but a 


ll 





small gradient must exist to prevent exhaust gas recircula- 
tion. Sufficient secondary air should be entrained to cool 
the engine exhaust to acceptable levels. Also, exhaust flow 
with very high kinetic energy can damage the test cell; 
secondary air entrainment reduces this energy. 

Augmentation ratio is defined as the ratio of entrained 
secondary air to the jet exhaust mass flow (primary flow). 
In practice, engine/augmenter tube separation is adjusted to 
generate sufficient secondary flow to prevent exhaust gas 
recirculation at all power settings. Occasionally, the 
augmenter inlet diameter is varied to further control the 
entrainment process. 

Both analytical and empirical methods have been applied 
to the entrainment mixing problem. Empirical equations are 
usually applied to particular test cells, but operator 
experience coupled with trial and error are required to 
create the desired secondary flow. ‘ 

Larger and more expensive test cells are needed in the 
future to replace obsolete cells and to accommodate new 
engines with greatly increased mass flow rates. Present 
design tools do not adequately describe the flow field in 
a test cell and new design concepts cannot be evaluated until 
the cells are constructed. An accurate flow model would 
provide a needed design tool and would help prevent costly 


design errors. 
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baer or ChiL POLLUTION 

Test cell installations are becoming increasingly sub- 
jected to localized state and community pollution control 
board regulations. The occasional prominent exhaust plume 
Ersing from test cells 1S a continual source of local com- 
plaint. Even when the plume is not visible, the cell emits 
pollutants similar to those of aircraft operating in the 
area. 

Fortunately, the stationary cell is amenable to exhaust 
treatment which has been shown to be effective in reducing 
particulate emissions to acceptable levels. The cost asso- 
ciliated with the construction and operation of the necessary 
equipment is high and is a function of total mass flow 
(primary plus secondary) through the cell. Thus, it is 
necessary to minimize total mass flow to minimize abatement 
costs. 

Carbon particulates and unburned hydrocarbons are of 
immediate concern in the abatement of emissions from jet 
test cells. Also of concern are nitrogen oxides (NOL), 
carbon monoxide (CO) and sulfur dioxide. The amount of each 
pollutant emitted is a function of engine type and mechanical 
condition. Power setting alters emissions with CO and unburned 
hydrocarbons predominate at idle. Particulates and nitrogen 
oxides are the main pollutants at high power settings. After- 
burner operation provides very high external exhaust tempera- 


tures and chemical reactions continue as the flow passes into 
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the augmenter tube. The secondary air mixing is believed 
to alter chemical reaction rates by effectively quenching 
portions of the exhaust jet. Present flow models do not 
describe the process completely and it has not been possi- 
ble until recently to accurately measure masS emission fac- 
tors in the afterburner core due to the high temperatures 
and the chemical kinetic nature of the process. 

To fully optimize the pollution abatement program it is 
necessary to understand how exhaust pollutants are formed, 
in what quantities and at what rates. The effect of engine 
power setting and secondary flow entrainment on emitted 
pollutants must also be analyzed. Pollutant meaSurement 
would be simpler and less costly if accurate species and 
velocity distribution predictions could be made at specified 
points in the test cell exhaust system. To accomplish these 
tasks, more complete knowledge of the air flow through test 


cells must be acquired. 


See oCOPE OF STUDY 

Most analysis for predicting dilution air is very empiri- 
cal, usually requiring an empirical jet spread equation. 
Recent advances have been made toward solving the recircu- 
lating, turbulent fluid flow problem. One technique 1s due 
Eomsoataimg et al. (Refs. 1 and 2]. The governing two- 
dimensional elliptic partial differential equations are 
reduced to finite difference equations and solved by numeric 


methods. 
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Air flow through a test cell is known to be turbulent 
and does possess regions of recirculation. The purpose of 
this study was to apply the methods described by Spalding 
to the flow field encompassing the engine exhaust and 
augmenter tube of a typical Navy jet engine test cell. 

Flow visualization and static pressure recovery were 
primary objectives of this study. Parametric studies were 
made which involved the variation of augmenter geometry, 
engine power setting, augmentation ratio and cell inlet 
flow conditions. 

The study was restricted to non-afterburning, subsonic 


engine exhaust conditions. 
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(ie eo RE at CAL APPROACH 


Two-dimensional fluid flow with regions of reverse-flow 
can be described mathematically by second-order elliptic 
partial differential equations. These equations can be 
reduced to "finite-difference" non-linear algebraic equations 
and solved simultaneously by an iterative procedure. One 
such procedure is detailed by Spalding and others in References 
1 and 2 and was the method used in this study. 

The conservation of mass, momentum and energy are the 
fundamental equations used in the problem formulation. The 
dependent variables chosen by Spalding et al. [1,2] are 
stream functions (wv) and vorticity (w) to assure mass and 
momentum conservation. Temperature was used as the dependent 
variable for energy conservation since the flow was considered 
inert. Since the velocities under consideration were low, 
kinetic heating was ignored. Specific heat was considered to 
be constant. 

Although pressure and, to a large extent, velocity are 
removed from the formulation, they can easily be "backed-out" 
of the vorticity and stream functions solutions respectively. 

A turbulence model is needed in order to provide 
“effective” transport properties. The Jones-Launder 
turbulence model [2,3] is used to describe the variation 
of the effective viscosity, Wore Two variables are used 


in this model: turbulence kinetic energy, K and turbulence 
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energy dissipation rate, €. The model is a time-averaged 
turbulence model and the two associated differential equations 
can be cast in elliptic form. Thus, the solution procedure 
for w, w and T also applies to K and e«. The relationship 


between K and ¢ is: 


Ke 


USet ne er 


where 


local density 


7 
ll 


empirically determined coefficient = 0.09 


or 
I 


These five equations together with the equation of state 
for a perfect gas describe the complete flow field. Choosing 
cylindrical coordinates (R,98,z) allows all five equations 
to be written in the following general form: 


(C4 >) 
g lag oxR! - 5ploqz)] - 5z PgR re eae 


5 a{c 0) 
where ¢ stands for Ve wir 6. The Lunctions Bor 
2a Cy and S 4 are given in Table I. 
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TABLE I 


COEFFICIENT VALUES 








db a b e S 
Ww 2 Z S 
i 
Y 0 1 a 
0 RZ R 
u 
K 1 ait 1 -(G - pe) 
K 
: , ae , 7 - Cope 
0 K K 
€ 
u 
T i 2st l 0 
h 
where 
U 
On = a taken to be 1.0 
K,eff 
Dy eff = the exchange coefficient of the kinetic 
: Snes soEmtiie cluctUating motiemeR: 
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oc, =———— taken to be 1.3 
a 7 
Get is 
r the exchange coefficient of the turbulence 
Gare tL energy dissipation rate, €. 
u 
On = —— taken to be 1.0 
h,eff 
er the exchange coefficient of heat 


oV dV V 
= Zee Ree eee 
ae et 2| (az cee aR” , 


ave OV, 9 
o (TR i OZ 


ea alG C5 are MiprnLecimyvecaeternned CoerLficients . 


a 


(c) = 125. C5 = Z0)) 


Neglecting tangential momentum, the vorticity source term 


ry 1s 
2 Z 
S -pf a (2 _ 83 
(w/R) oz 2 oR 
0 v," a Va op 
ae Soa) 32 | Tene dilcg sole terns 
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The first term of the general equation is the convective 
term while the second and third terms are the diffusion terms. 
The last term is the source term which includes all terms not 
strictly convective or diffusive. 

Once the equations are solved and a converged solution 
is acquired, the static pressure distribution in the flow 
field can be calculated by integrating the momentum equation 
in the axial or radial direction. The resultant equation for 


the gradient in the axial (z) direction is 


" QV 
dP _ 9 a ee Z 2 
2 d2 ) eff E oZ 3 div 7 — | 


5 ove OV, 
OR Seleccsl sms ° oe 7 ReVIYR 


The analogous form in the radial (R) direction (neglecting 


+ 


Zh 


@-momentum) is 


2V 
2 ae cee 
aR E | voce (= 2 aiv @) 


where 
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Pressure is then obtained by the following simple 
integration formula (which assumes a linear variation in 
pressure gradient between grid points): 

ii oP dz nie 1 F(2e + ee ] (2 - 2 ) 
OZ OZ 


a OZ 2 A B B A 


Zi 





ESE. MODEL DISCUSSION 


Meee COORDINATE SELECTION 

Although the test cell cross-section is baSically rectan- 
gular, the engine exit and augmenter tube can be accurately 
described using two-dimensional axi-symmetrical flow. Of 
primary interest is the mixing of the engine exhaust (primary 
flow) and ambient air (secondary flow). This process pri- 
Marily occurs between the two circular "tubes" and within 
the augmenter tube. Reference 4 suggests that test cell 
wall proximity will influence primary/secondary flow mixing. 
However, the large dimensions of the cell compared with the 
engine and augmenter diameters should allow the use of 


cylindrical coordinates with little error. 


Eee) GEOMETRY 

Figure lis a profile view of a test cell located at 
NARF Alameda, used to test TF-4l1 engines. Figure 2 is the 
Pere try established for an engine test while Figure 3 is the 
geometry used in this study. The actual dimensions are 
easily varied in the computer program allowing parametric 
studies to be made to assess the effects of geometry change 
on the flow field. The field 1S assumed axi-symmetric although 
in practice the engine/augmenter combination is located 
nearer the test cell floor than the walls or ceiling. The 
cell ceiling is modeled as cylindrical but in actuality 1s 


amate flat. 
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TO minimize the number of grid points required in the 
axlal (z) direction, the cell "inlet" was located a few feet 
upwind of the engine exhaust. 

The straight augmenter duct "lip" used in the computer 
model was a necessary simplification to the beveled lip 
Since grid points in the coordinate system lie along mutually 
perpendicular grid lines. Augmenter tubes in Navy test cells 
are by no means uniform in design and inlet shapes vary 
greatly. It is believed that a straight lip design does 
not greatly affect the flow analysis. 

Augmenter tubes used to test afterburner operation 
usually contain quench systems and flow baffles. These were 


not considered in this analysis. 


C. GRID SPACING 

Various grid sizes were used to appraise the effect of 
grid size on the predicted flow field. The largest number 
of grid points was described by a 60x50 (zxR) matrix and te 
smallest by a 20x25 matrix. Regardless of the size of the 
grid, every boundary must be defined by a grid line. [In 
general, grid line separation should be smaller as one 
approaches a wall since gradients in properties are large 
in this region. 

Spalding [1,2] suggests that in the region of solid 
boundaries each succeeding grid spacing not differ from the 


Hast by more than a factor of 1.5. 
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D. BOUNDARY CONDITIONS 
All dependent variables must be specified in some manner 
on the boundaries. The boundaries of this model include: 
(1) the engine exhaust duct, (2) the cell inlet, (3) the 
solid boundaries, (4) the augmenter outlet, and (5) the 
geometric centerline. 
iPeencine Exhaust Duct 
a. Stream Function/Velocity 
Plug flow was assumed at the engine exhaust 
plane. The axial velocity distribution is uniform (Figure 
4) to produce a specified mass flow and no radial component 
Srnvelocity exists. 
Dremee vO Eleisty 
Since no velocity gradients exist, vorticity 
is identically zero. 
c. Mass Fraction/Temperature 
A non-dimensional mass fraction/temperature was ' 
assigned a value of 1.0 at the engine exit. This corresponds 
to the actual exhaust temperature existing at the specified 
power setting. 
ad. Turbulence Kinetic Energy 


Kinetic energy of turbulent motion (K) is defined 


by 
vp tv) + vy 


Wie Gey are the instantaneous fluctuations 
in free stream velocity. 
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For isotropic turbulence the expression for K 


Cammoe written K = 5 a 


Intensity of turbulence is defined as 


a 2 2 2 
z (Vp oo Vo + Vv, ) 


<I 


Assuming isotropic turbulence 


where V is the mean flow velocity. 


Spalding [1,2] presents an empirical relation 


herek 1f Uniform flow: 
ae 
K = 0.004 (V) 
Setting 


a =e UIA co 


bo} We 


and substituting into the intensity relation we find 


-008 =2 
3 V 


V 
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schlichting [5] experimentally determined 
intensity values of between 0.05 and 0.08 in low subsonic 
wind tunnel flow. Thus, it appears Spalding's relation is 
adequate for estimating K in the engine exit and was utilized 
as the boundary condition. 
e. Turbulence Energy Dissipation Rate 
K and € are related through the dissipation 


Hength scale £ [2] by 





In this program £ was taken to be: 


£ = 0.022 x (Engine Radius) 


aeotakec to Cell 
a. Stream Function/Velocity ‘ 
A linear variation of inlet axial velocity was 

chosen which varied from a positive value at the engine 
wall to zero at the cell ceiling to produce a specified 
secondary mass flow (Figure 4). No specific information was 
available to specify this flow distribution accurately. 
Lacking more specific information, the radial components of 
velocity (V,) were fixed at zero. 

DreeevOEETCILEY 


Vorticity was computed using 


AE 
AI 





c. Mass-Fraction/Temperature 
Non-dimensional temperature and mass-fraction 
were asSSigned a value of 0.0 at the cell inlet. fThis 
corresponded to an antual temperature of 520 °R., 
dad. Turbulent Kinetic Energy (K) 
K was specified at each grid point using the 


expression 


a0 00 4, 


Where Ta is the local flow velocity. 


e. Turbulence Dissipation 


gE = K>/2 72 


cs 
I 


0.022 x (Cell Radius - Engine Radius) 


3. Solid Boundaries 

ti wemareseignt seid bounGarres Gescribed in the 
present model. Each is numbered in Figure 5 for reference. 
Two are engine surfaces, two are cell walls and four describe 
the augmenter tube. 

a. Velocity/Stream Function 

Along the engine surfaces (1 and 2, Figure 5) 

stream function is constant and determined from the specified 


mass flow through the engine. The remaining solid boundaries 


a 





have the stream function fixed at a value which is the sum 
of the engine mass flow and cell inlet mass flow. All 
velocities at a solid boundary are zero. 

Dae VOGtLer ty 

Vorticity near a solid boundary was computed in 
one of two ways depending on the magnitude of the local grid 
cell Reynolds number. 

(1) Laminar Sublayer Case. If the local grid 
cell Reynolds number (R = Ay/uR) is found to be less than 
290 [1,2], a linear velocity gradient is assumed. Using the 
additional assumption that the velocity gradient in the 
direction perpendicular to the wall is much greater than the 
gradient parallel to the wall, the vorticity is calculated 
uSing the vorticity definition and discarding the parallel 
term. For example, given a wall parallel to the horizontal 
(z) axis with flow in the positive z-direction, the vorticity 


would be found by: 





7 ee ae ee uf 
GU) eee ee Se 
OZ oR oR AR AR 
oV 
where xD is considered small. Vp is the velocity in the 
z-direction at node P (Figure 6). Since velocity at the wall 


is zero, it is easily seen that Vp is twice the average 
velocity between the wall and the nearest grid point off the 
wall. This is convenient since the average velocity is 


related to the stream function by 
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Therefore, the vorticity expression is given by 


Cc 
< 
< 


- a ee os RC 
= lla atl reee 


o (AR) 


which is easily calculated. 
(2) Turbulent Case, When the local grid cell 
Reynolds number at the near wall node exceeds 290 the flow is 


x 


considered turbulent. Vorticity is then defined as 


"shear 
dW) = a eS 
Hwall 
where ay lis the wall shear stress and Wall is the 


Viscosity on the wall. 1S given by the relation: 


“shear 


_ Sf 92 
Tshear 92 ? 


c 
The coefficient (>) was calculated using 


cs eg. ey fo ae 
2 Re'*> u Ry 


c 
ae where Re = a 
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(3) Special Points. The preceding discussion 
of vorticity calculation near solid boundaries assumed a 
well-developed boundary layer. There are locations in the 
field where the fluid may meet a parallel solid surface and 
the vorticity calculation becomes suspect. Such a point is 
the intersection of surface 5 and surface 6 of Figure 5 
(point C). THere are two problems associated with this point. 
If the flow is predominantly in the positive z-direction, a 
fully developed boundary layer will not exist on surface 5 at 
point Cc. If the flow is in the negative R-direction near 
surface 6, the identical problem arises. Also, each surface 
(5 and 6) will have a unique value of vorticity at point C 
and it is necessary to keep them properly identified in the 


computer program (see sketch below). 


surface 5 


6 
D [ 
Z 


No completely satisfactory numerical method 
exists in this scheme to account for the singularity. Initial 
computer runs were made allowing the vorticity at point C 
to be calculated as if a fully developed boundary layer 
existed on both surfaces (5 and 6). Converged solutions 


resulted with realistic flow visualizations around this point 
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as long as a recirculation region did not develop on the 
eugmencer Guct lip (surface 6) . 

When flow developed in the positive R-direc- 
tion near surface 6 with negative z-direction flow on surface 
9, the solution would not converge and often diverged (see 


sketch below). 


surface 5 


@ 6 @ 


Surface 
6 





References 1 and 2 state that special 
attention must be paid to points which define enlargements 
in pipe flow (Points A and D of Figure 5). In order to 
obtain realistic jet spread, stream function must be fixed — 
on the node immediately downwind of these points. It became 
evident that the flow would usually be in the negative z- 
Girection along surface 5, so wy was fixed at the node circled 
in the sketch above. This solved the divergence problem 
occurring under certain flow conditions, but appeared not to 
aterm the solution appreciably. Point ® is unusual in that 
it defines the end of an “infinitely thin" wall. It also 


defines an "enlargement" and for this reason, the immediate 


Sal 





downward WY value was fixed at the value aSsigned to 





point D. 
surface 8 
Bre aice surface 
6 i 
@ ele e 
oe ele el fixed 
D 
e o € 
F 


dominant flow 
direction 


For purposes of calculating the vorticity at node F (see 
above sketch), the vorticity at point D was taken to be the 
average of the vorticities computed on surface 6 and 7 shown 
in the above sketch. 

c. Turbulence Kinetic Energy 

Turbulent flow produces steep gradients of all 

transport properties near walls. To reduce the required 
number of grid points and to reduce computation time a 
"Slip value" of turbulence Kinetic energy (K) is calculated 
on the wall. The “slip value" is fixed to produce an arith- 
metic-mean gradient between the wall and adjacent node, which 
equals the gradient of the Couette flow solution midway 
between the nodes. A complete analysis of this method can 


be found in References 1 and 6. 
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ad. Turbulence Dissipation 
The effective boundary for turbulence dissipa- 
tion (€) is taken to be the n grid points adjacent to the 


wall. Wolfshtein [6] presents the value at the near-wall 


node as: 
“u 3/2 
a ee) ase 
u 
where Ep = near-wall value 

7 = empirical constant = 0.09 

a = empirical constant = 0.22 

Kp = near-wall turbulence kinetic 
energy value 

AR = wall to near-wall node 
distance 

e. Temperature/Mass Fraction 


All solid boundaries were assumed adiabatic 
and, of course, impermeable. Therefore, no gradients in 
temperature or mass fraction existed between the walls and 
near-wall nodes. 
4. Augmenter Exit 
Axial gradients of all dependent variables were 
taken aS zero. 
5. Geometric Centerline 
Radial gradients of w/R, K, e€, T and mass fraction 
are taken as zero. Also the centerline is the stream function 


fim O) reference line. 
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E. INITIAL CONDITIONS 

The flow field must have initial values assigned to all 
variables. Experience is required to fill the field with 
compatible variables to prevent divergence in the first few 
iterations. 

It was found that a reasonably accurate stream function 
distribution was essential to prevent divergence. In par- 
ticular, the predominant flow direction must be established 
and large discontinuities in stream function should not 
exist. The flows in the augmenter tube and the free stream 
area between the engine exhaust and augmenter tube were 
initialized as well-developed pipe flow. This allowed com- 
patible vorticity values to be calculated in these regions. 

Density, mass fraction, viscosity, and turbulence 
kinetic energy were specified using reference values taken 
from either the cell inlet or engine exit. Turbulence 
dissipation was calculated using the relationship discussed 


under boundary conditions. 


Peo OLULTLON PROCEDURES 
To counter computational instabilities in the first 400 
iterations, stream function variation between iterations was 
under-relaxed by a factor of 0.3. All other variables were 
under-relaxed by 0.5. The relaxation parameters were 
increased to 0.6 and 0.8, respectively, after 400 iterations. 
Various flows, geometries, and grids were investigated 


successfully using the above relaxation parameters. They 
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precluded divergence in the most unstable flow condition 
handled: high engine mass flow with low augmentation ratio 
(200 lbm/sec, 0.25, respectively). However, it should be 
motea that initial conditions must be well established for 


these relaxation parameters to be effective. 
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IV. MODEL DEVELOPMENT 


The geometry shown in Figure 3a was chosen to represent 
a test cell section containing an engine exhaust nozzle and 
augmenter tube. Approximate dimensions were taken froma 


TF-41 turbofan engine undergoing testing at NARF Alameda. 


A. GRID SPACING 

A grid composed of twenty lines in the z-direction and 
twenty-five lines in the R-direction (20x25) was chosen to 
initially describe the flow field. (Intersections of 
eewainy perpendicular grid lines define grid points.) 

This number of grid lines allowed the use of small grid 
spacing near solid walls but produced relatively large 
spacing away from walls. 

Computation time is directly proportional to the number 
Pmumcianemperntes in the field. A small number of grid lines 
was chosen to minimize computer time during program FeRtccrtian 
Three hundred iterations typically required 7.5 minutes 
(excluding compile time). It should be noted that vorticity 
source terms were ignored during this initial period. 

Vorticity source terms and pressure recovery calculations 
were added only after the program was successfully iterating 
to a converged solution. The number of grid lines was 
inereased to 40x40 (zxR) to allow for pressure gradient 
calculations without the use of wall or near-wall variable 
values. OtherwiSe, large errors accumulated in the pressure 


integration process. 


36 





Eleven hundred iterations typically required 109 minutes 
of computation time (excluding compile time). This time 
includes one pressure calculation along eight lines in the 
flow field and those computations performed in the flow 
initialization process. However, the total time for the 


latter calculations amounts to less than one minute. 


B. CONVERGENCE CRITERIA 

Convergence rate was influenced by flow conditions, grid 
spacing and relaxation parameters. Two criteria were used 
to appraise the convergence level of the solution. 

Piety eene Percentage ciange Ing, Wp, K, 6 and mass 
Piaadet ron (at all grid points in the field) must be arbi- 
trarily small between succeSsive iterations. This was done 
computationally by storing the largest change occurring in 
the field for each variable and printing that value after 
each iteration. 

The second convergence criterion proved to be the most ° 
conclusive. That is, when pressure calculations along each 


grid line remained unchanged (with increased iterations) 


convergence was assumed. 


C. PLOTTING PACKAGE 

Solution values of stream function and temperature were 
transferred to computer cards at the conclusion of each run. 
These cards were then used as input to a second program 
which plotted streamline contours and lines of constant 
Eemperaeure On Separate graphs. A CALCOMP Model 765 plotter 


Was utilized for this purpose. 
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Veer at Cre TUDY 


Geometry and inlet conditions were varied to appraise 
model utility and sensitivity to flow conditions. TF-41 
engine test data were used to establish engine exhaust static 
temperature and velocity. Static pressure at the engine 
exit was fixed at atmospheric pressure (2116.0 lb -/£t*) for 
this study. Engine exit velocity was found by computing p 
from the equation of state for a perfect gas and using the 


one-dimensional mass conservation equation. 


and 


™m . ; 
engine 


V 
engine 0 Aangine 


where m total mass flow through engine 


(including bypass air) 


A = area of exhaust plane 


Engine exit static temperature was computed using the 
turbine outlet temperature and the bypass ratio given in 
Reference 7. A summary of these parameters is shown in 


Table II. 
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TABLES. tT 


te=4)] Test CONDITIONS 


——_—_—— LTS CC NT CC A fA 





nc m (lb /sec) De (oR) Voxit (ft/sec) Boa 
Number 

Idle 100 800 594 0.43 
802 LPS6 900 

90% 200 927 W376 0292 
94% 238 955 
97% 249 970 
100% 263 1001 


A specific flow condition was chosen to provide a solution 
to which all succeeding solutions could be compared. That 
condition was: idle engine power setting with an augmentation 
ratio (AR) of 0.5. The geometry was fixed as shown in 
Figure 3. 

The "reference" flow condition was chosen as it provided 
low subsonic flow conditions at the engine exhaust plane 
(M = 0.43) and the augmentation ratio was typical for actual 


mest COnGLtions. 


Bo 





Pressure calculations were made along the dashed lines 
shown in Figure 7. The lines show only approximate locations 
Since grid spacing was often altered when geometry was 
varied. 

Eight computer runs were made to examine the usefulness 
of the model. Each run was allowed approximately 2.5 hours 
(1400 iterations) to assure solution convergence although 
most runs achieved practical convergence prior to this 


time. A summary of the study follows in Table III. 


TABLE III 
FLOW CONDITIONS 
ke Reference flow condition was established 


as discussed above: Idle Power, AR = 0.5, 
Figure 3a geometry. 


oe Two runs were made varying the AR: fIdle 
Power, AR = 0.75 and 1.0, Figure 3a 
CeOmMe tive. 

oe Two runs were made to change the length 


of the augmenter duct lip: Idle Power, 
AR = 0.5, Figure 3b, geometries. 


4, One run varied the engine to augmenter 
separation distance: Idle Power, AR = 0.5, 
Figure 3c, geometry. 


or. One run varied the augmenter radius: 
Idle Power, AR = 0.5, Figure 3d, geometry. 


6. One run changed the power setting of the 
engine: 90% power, AR = 1.0, Figure 3a 
geometry. 
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Vie Senco lb wom ly Dl oeuUSo LON 


A. FLOW CONDITION 1 (Reference Condition) 

Streamlines and static temperature contours are shown 
in Figures 8a and 8b, respectively. Pressure gradients along 
axial and radial grid lines (Figure 7) are graphically 
displayed in Figures 8c through 8f. Pressure contours along 
aXial and radial grid lines are displayed in Figures 8g 
through 81i. Figure 8} presents the variation in augmenter 
pressure rise as a function of iteration number. 

1. Streamlines/Temperature Contours 

Two-dimensional flow solutions are produced by this 

study. Therefore, the flow field consists of a one-radian 
"portion" of the axisymmetric field (Figure 8a) and the 
streamlines are labeled accordingly. For example, the mass 
flow rate of a TF-41 engine is 100 1b, /sec at idle power 


setting. The two-dimensional mass flow would be 


its - Moctual 
2-D 27 


100 lb_,/sec 
se eee ee S59 lb_/sec-radian 


271 
aos shown ain Figure sa. 
Two distinct regions of recirculation are apparent 
in Figure 8a. One prominent region exists above the augmenter 
tube bounded on two sides by the test cell ceiling and wall. 
Points of stagnation are shown located on the exterior 


augmenter tube surface and cell wall. 


Al 











The second prominent region of recirculation is 
located downwind of the augmenter inlet lip. Close exam- 
iiaeren Of the streamline having its origin at the lip 
reveals axial flow for a short distance followed by rapid 
directional change. Axial flow could have been the result 
of fixing the stream function at the first downwind corner 
as discussed in Section III.D.3.b.(3). However, an additional 
computer run was made with the same flow conditions except 
that the stream function at the downwind node was allowed to 
be computed normally. The solution for stream function was 
essentially identical to the first run. In any case, the 
lip has only small effects on the shape and size of the 
recirculation region. The recirculation region is formed 
at the low thrust, low augmentation ratio condition pri- 
marily by the radial component of the augmentation air as 
it flows into the augmenter. 

The region of flow reattachment is well defined on 
the augmenter interior surface and is located approximately 
two augmenter diameters downstream of the inlet. 

Temperature contours shown in Figure 8b are an 
indication of the amount of mixing occurring between 
engine flow and secondary air. From the non-dimensional 
mass fraction and temperature definitions, a thoroughly 
mixed flow would have a static temperature of 707 °R. 

Under these flow conditions, a well-mixed flow occurs 
approximately four augmenter diameters downstream 


of the augmenter inlet. The temperature contours 
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clearly show the effect of secondary air entrainment and 
recirculation within the augmenter tube on the temperature 
field. 
eee Pressure Recovery 
Sti eekly (vatelabialexeisalioy ql), 

Pressure gradients along lines A-A through E-E 
(Figure 7) are shown in Figures 8c, d and e. Lines labeled 
A-A and B-B are axial lines having their origin at the 
engine exit plane and terminating one node upstream of the 
augmenter exit. Lines C-C, D-D and E-E extend between the 
third node in the axial direction and the node immediately 
upstream of the augmenter exit. 

Large negative pressure gradients exist along all 
lines from the engine exit plane to just inside the augmenter 
inlet plane. Lines A-A and B-B exhibit the largest gradients 
in this region and examination of the solution reveals the 
exhaust jet rapidly departing from the uniform inlet velocity 
(plug flow) condition. Secondary air 1s being entrained in 
this area with line B-B being nearest the turbulent mixing 
zone of the exhaust jet. This may account for the largest 
pressure gradients existing along line B-B. 

Pressure gradients along axial line F-F (Figure 7) 
are very small and are not shown. Static pressures along 
axial lines A-A through E-E are shown in Figures 8h and i. 


These pressure profiles are all initialized at 2116.0 


2 
lb, /ft : 
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b. Radial (R-direction) 

Pressure gradients and static pressures along 
G-F are shown in Figures 8f and 8g, respectively. The engine 
exhaust jet turbulent mixing zone is clearly indicated by 
Fluctuating pressure gradients (Figure 8f). Static pressure 
is seen to increase toward the engine exhaust jet centerline 
which is due primarily to the secondary air momentum in 
the negative R-direction. 

No appreciable pressure gradients existed along 
radial line H-H which is located downstream of the flow 
reattachment point. 

c. Total Axial Static Pressure Rise 

Total rise (AP) was defined as the static pressure 
difference between end-points of the lines discussed above. 
AP is plotted against number of iterations in Figure 8j. The 
upwind end pressure for all lines is reference pressure but 
each line is seen to exhibit a unique value for AP at 1400 
iterations. No pressure gradients exist in the radial 
direction downwind of the reattachment point. Therefore, 

AP values should be identical in magnitude for lines A-A 
through E-E when the solution is converged. 

Close examination of the pressure gradients 
(Figures 8c through 8e) reveals line E-E to have small, "well 
behaved" gradients compared with lines A-A through D-D. 
Figure 8j indicates that the AP calculation for line E-E 
"converged" relatively early (1000 iterations). For these 
reasons, pressure rise was taken to be that computed for 


E-E (47 lb/ft’). 
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The differences in the pressure rise calculated 
along axial grid lines can be due to several things. The 
number of iterations has a dominant effect on the calculated 
pressure rise. The dependent variables (w, yw, K, €, T) vary 
only slightly with increased iteration time. However, the 
pressure calculation is apparently very sensitive to small 
changes in the stream function and density. This is probably 
because ab is a function of second derivatives of wv. In 
addition, pressure rise is calculated by assuming linear 
Variations in gradient between grid points and pressure has 
been held fixed for purposes of calculating the density 
(1.e€., 9p = p(T)). Improved accuracy might be obtained if 
pressure were calculated periodically in the solution in 
order that pressure effects on density may be taken into 
@ee@eunt. Lt was found that finer grid spacing did not 
appreciably affect the pressure calculation. Another possible 
error source may enter through the use of upwind differences 
and/or the inadequacy of the Jones-Launder model for the flee 
geometry under investigation. However, checks on the 
solution indicated that mass was conserved at each axial 


location with less than 14 percent error. 


B. FLOW CONDITION 2 

Augmentation ratio (AR) was varied with idle power setting 
forwassess the effect on total static pressure rise through the 
Eyoecem. Solution convergence and pressure gradients were 


famidar to Flow Condition 1 and axial line E-E was again 


45 





chosen to determine the value of AP. The solutions along 
the other axial lines had approximately the same slope on an 
AR vs. AP plot but were different in magnitude as discussed 
above. 

Figure 9a shows decreasing AP with increased AR. This 
result qualitatively correlates with Bailey [8,9] who experi- 
mentally varied exit pressure on a small jet pump and deter- 
mined secondary mass entrainment. As the augmentation air 
is increased for a fixed engine flow there is less momentum 
differential between the two streams. This apparently causes 
the decreased pressure rise with increased AR. In actual 
test cell operation there 1S a unique AR for each engine 
flow rate (assuming the augmenter is fixed in design and 
position). The actual value is determined from the pressure 
rise in the augmenter tube being equal to the pressure losses 
through the rest of the test cell. If pressure losses 
through the inlet and exit stacks (with eyed associated 
acoustic and pollution abatement equipment) and through the 
horizontal portion of the cell can be estimated as a function 
of flow velocity, then the AR can be found in an iterative 
Manner by using the AP vs. AR plot determined with the model 
for a particular augmenter design and position. 

Figures 9b and 9c show the effect of AR on the stream 
function and temperature distributions, respectively. As 
the AR is increased the axial velocity component of the 


secondary air dominates over the radial component. The latter 
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component was observed above to be the major cause of the 
recirculation region. Thus, as the AR is increased, the 
recirculation region within the augmenter is seen to decrease 


in size and mixing occurs more rapidly. 


ce. FLOW CONDITION 3 
The augmenter duct lip length was varied by 0.63 ft 

(2.1 ft to 2.73 ft augmenter inlet radius). Stream function 
and temperature fields are shown in Figures 10a through 10d. 
No appreciable change in calculated pressure rise was observed. 
This result apparently occurs because at low augmentation 

ratios the size of the inlet lip has very little effect on 
the size of the recirculation region within the augmenter. 
A greater effect of duct lip length on pressure rise could 


be expected at higher augmentation ratios. 


D. FLOW CONDITION 4 

Engine to augmenter inlet separation was increased to 
ost (trom 1.5 £t). “Stream function and temperature 
contours are included as Figures lla and 1llb. No appreciable 
change in pressure rise was noted along axial grid line E-E 
when compared with Flow Condition l. 

The stream function solution indicated that no appreciable 
change occurred within the augmenter when the distance was 
imeredsea. Again, a more Significant effect could be expected 
for higher augmentation ratios. Moving the augmenter closer 


to the engine should have a major effect on the recirculation 
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region and the calculated pressure rise. As discussed above, 
the inlet velocity profile for the augmentation air was 
assumed to increase linearly from the cell ceiling to the 
engine at a location just upstream of the engine exhaust. 
This profile would change as the augmenter position is changed. 
Since test cells are designed to provide uniform flow at the 
engine inlet, future investigations with the model should be 
Maceewith the secondary air inlet specified as uniform at the 
engine inlet plane. In this case, the velocity profile near 
the engine exhaust plane would be free to develop in a 


meattstic Manner. 


Ewe ELOW CONDITION 5 

Augmenter radius waS increased by 0.5 ft to 3.5 ft while 
keeping the duct inlet diameter fixed. Figures 12a and 
12b present the temperature and stream function contours. 
Pressure rise (AP) along E-E was diminished to 37 yey ies 
from the reference condition value of 47 eect Although 
only one run was made, this pressure loss is considerable an 
indicates AP decreases with increasing augmenter diameter 
(given fixed inlet flow conditions). The increased size of 


the recirculation region (which consists primarily of 


augmentation air) May contribute to this pressure loss. 


F. FLOW CONDITION 6 
Engine power setting was increased to 90% and the A R 
was set at 1.0 which is typical for an actual engine test. 


Figure 13a shows the recirculation region smaller with flow 
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reattachment farther upstream. The temperature contours 
(Figure 13b) show engine exhaust and secondary air being 
thoroughly mixed farther down the augmenter tube than 
occurred at idle flow. 

Pressure rise versus number of iterations is shown in 
Figure 13c. The numerical values of the dependent variables 
did not change appreciably after 1400 iterations. However, 
AP's differed in value considerably along lines A-A through 
E-E, so no definitive pressure rise could be derived from the 
results. The solution procedure (which assumes incompressible 
flow and uses the associated elliptic equations) used in 
this study is questionable when applied at the higher subsonic 


Mach numbers. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


ee CONCLUSIONS 

The model provides two-dimensional flow visualizations 
in turbojet test cells and augmenter tubes for low subsonic 
flow conditions. These can be used to identify regions of 
recirculation and to assess the amount of mixing occurring 
between engine exhaust gases and secondary air. Optimum 
locations for pollution sampling equipment can be selected 
by examining the numerical solutions. 

Static pressure rise within the augmenter tube can be 
calculated for a given flow condition and augmenter design 
from converged solutions, provided that regions of "steep" 
gradients in variables are avoided. Total pressure rise 
through the system can be found and factors which influence 
this rise (augmenter design and location, engine flow rate, 
test cell size, and etc.) can be studied. 

High subsonic flows invalidate the use of partial 
differential equations of the elliptic type; therefore, 
the model provides realistic solutions only at low subsonic 


Mach numbers. 


B. RECOMMENDATIONS 

1. Additional computer runs should be made to firmly 
establish the effects of augmenter and cell geometry on flow 
mixing and pressure rise (AP). For example, the effects of 


moving the augmenter inlet closer to the engine exit plane 
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and moving the test cell wall forward to the augmenter inlet 
plane should be considered. In addition, the program could 
also be easily modified to extend the augmenter duct lip in 


the opposite direction. 


2. The cell "inlet" position should be moved upstream 
to the engine inlet plane by the addition of a few grid 
lines in the z-direction. Plug flow could be assumed at 
the "new" cell inlet and the resultant flow field would show 
the actual velocity distribution at the "cell inlet" used 
tempthis Study. in addition, this would permit calculation 


of exhaust gas recirculation. 


3. Static pressure and velocity measurements should be 
taken under various test conditions at NARF Alameda. The 
Size and shape of the recirculation region located above 
the augmenter should be qualitatively verified by photo- 
graphing a tufted grid placed between the test cell augmenter 
tube and ceiling. These data would help in validation/ 


improvement of the computer model. 


ou 





ASLENOSD 113) LSal 





40399] 4ap 


9qn}z uezuewbne 


dLasn09e 
qsneuxa 


[| [99 


“L sYyndI4 


autbua 


SL aued 
2137Sn05e 


ae 
inl 


yaLuL 
[a0 


52 





NOILVYNDTANOD LS3l Lo-sl 


pono WY Bl 


a 
wo 


Pe ee ee Gt 
ee 4a2uawbne 





4 


[Lem 
L129 





“¢ dang 








SS 


W/W 


auLbue 





a a oleae eae oe 


v Ciao 


oo 


But} load Le 


EL 


353 





T NOILIGNOD MOT 
>AULSWOS9 TSCOW YSLNdWOD “PE JUNIIS 


Pa $e 0 
1"9- —_,— - 


ysneuxea . 
| auLbua ns el . 


——— 
| SOP y ais auLbua 
A3zuawbne 

| 


S wOtby FONP 
a, Aazuawbine 


adeJuns daqzuowbne 








gr en ee NT TTT mee 





aa 


ee ~~  e me 





OL | 

LLeM | 
Bg. OL (sasodund 

[apow 40s) | 

FaLUl | 
LLao 

| 

| 


But, Lad | 19a 


54 





| augmenter | augmenter 


eS Sr 








ee] , OSI0), 
3 t 
| engine ngine 
a) e. = 
sce is t<— 1.5! | 
C.L.-———- - —_—_— - ———-C.L. 
large augmenter inlet sma}] augmenter 
radius inlet radius 


FIGURE 3b. COMPUTER MODEL GEOMETRY: FLOW CONDITION 3 
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FIGURE 3d. COMPUTER MODEL GEOMETRY: FLOW CONDITION 5 
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FIGURE 5. IDENTIFICATION OF BOUNDARIES 





solid boundar 
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FIGURE 6. LAMINAR SUBLAYER VELOCITY PROFILE 
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FIGURE 8f. RADIAL PRESSURE GRADIENTS: 
FLOW CONDITION 1 
(LINE G-H) 
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FIGURE 8g. RADIAL STATIC PRESSURE VARIATIONS: 
FLOW CONDITION 1 
(LINE G-H) 
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FIGURE 8j. AP RISE VERSUS NUMBER OF ITERATIONS: 
FLOW CONDITION 1 
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FIGURE 9a. AP RISE VERSUS AUGMENTATION RATIO: 
FLOW CONDITION 1 
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FIGURE 13c. AP RISE VERSUS NUMBER OF ITERATIONS: 
FLOW CONDITION 6 
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